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Several emergent phenomena and phases in solids arise from configurations of the electronic 
Berry phase in momentum space that are similar to gauge field configurations in real space such as 
magnetic monopoles. We show that the momentum-space analogue of the “axion electrodynamics” 
term E • B plays a fundamental role in a unified theory of Berry-phase contributions to optical 
gyrotropy in time-reversal invariant materials and the chiral magnetic effect. The Berry-phase 
mechanism predicts that the rotatory power along the optic axes of a crystal must sum to zero, a 
constraint beyond that stipulated by point group symmetry, but observed to high accuracy in classic 
experimental observations on o-quartz. Furthermore, the Berry mechanism provides a microscopic 
basis for the surface conductance at the interface between gyrotropic and nongyrotropic media. 
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The topological consequences of time-reversal symme¬ 
try breaking in two-dimensional electronic systems have 
been a focus of interest since the discovery of the quan¬ 
tum Hall effects [T]. Similarly interesting phenomena 
arise from breaking inversion symmetry (IS) in three- 
dimensional systems, for example in one type of ”Weyl 
semimetal, ” 013], possibly realized in TaAs 00, where 
IS breaking allows for non-trivial topological states that 
contain pairs of chiral gapless fermions. At least in insu¬ 
lators, it is now widely known that there exist quantized 
transport phenomena as a result of topological invari¬ 
ance. One goal of this paper is to demonstrate an exam¬ 
ple of topology in the optical response of metals , which 
can be derived using relatively simple semiclassical elec¬ 
tron motion in low-symmetry solids. 

The main effect to be discussed, natural optical ac¬ 
tivity, arises in materials that break IS but retain time- 
reversal symmetry. We find two unexpected features of 
optical activity, one of which may already have been ob¬ 
served in <a-quartz, and obtain a general constraint on 
optical activity in some frequency ranges that impacts 
some designs for topological photonic devices. Electron 
dynamics in such materials is subtle because, despite the 
lowering of spatial symmetry, the energy spectrum it¬ 
self remains symmetric, z.e., e(k) = e(—k). Thus, the 
physics that underlies transport anomalies in such sys¬ 
tems must involve the properties of the electronic wave- 
functions themselves, rather than their energy levels. 

It is by now generally understood that the 
wavefunction-dependent transport properties of electrons 
on a lattice are affected by the Berry curvature, ffc(k), of 
the Bloch states [MB- In the presence of a nonzero 
fi(k), the semiclassical equations of motion for an elec¬ 
tron wavepacket are modified to respect the duality be¬ 


tween position and momentum space, 

r(k) = v(k) + kx Cl(k) 

—{h/ e)k(r) = E(r) + i x B(r), (1) 

where v(k) = /i _1 V k^(k), Cl = V x (uk | — and 

electron charge is (—e). From the symmetry of the mod¬ 
ified equations of motion with position and momenta, it 
is clear that Cl (k) can be viewed as an effective magnetic 
field in momentum space. In this paper we introduce an¬ 
other useful momentum-position space correspondence, 
involving the dual to the magnetoelectric coupling term 
in the Lagrangian density, that is, Cij{r,t) = otijEiBj. 

The dual nature of the semiclassical equations suggests 
a corresponding tensor in momentum space, 

0«(k) = i>i(k)nj(k), (2) 

which turns out to play a fundamental role in a uni¬ 
fied theory of Berry-phase contributions to the transport 
and optical properties of inversion-breaking media. The 
scalar diagonal part of is topological and referred to 
as “axion electrodynamics” mm, and the trace of Gij 
also has a topological significance in multiple contexts. 

As is clear from Eq. [ljthe signature of nonzero Cl (k) 
is the existence of an ” anomalous” current transverse to 
the applied force. For ac electric fields the transverse cur¬ 
rent manifests as the phenomenon of optical gyrotropy, 
in which a medium exhibits a different index of refrac¬ 
tion for left and right circularly polarized light 05]. Gy- 
rotropy in media with broken time-reversal symmetry 
is known as Faraday rotation, whereas in time-reversal 
symmetric systems it is usually referred to as ” natural 
optical activity” (NOA). Below, we develop a semiclassi¬ 
cal theory of NO A originating from the Berry curvature 
of inversion-breaking media, obtaining two new results. 
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First, the gyrotropic tensor, < 7 ^, that emerges from the 
topology of the Berry curvature is traceless. This rep¬ 
resents a constraint on the components of gij beyond 
well-known relations imposed by the point-group sym¬ 
metry and is therefore a signature of the Berry phase 
mechanism. Tracelessness of the topological poten¬ 
tially resolves an 80 year old mystery concerning NO A in 
<a-quartz and other materials on on. The second result 
is the existence of a surface current that flows in response 
to an electromagnetic wave incident at the interface be¬ 
tween gyrotropically active and inactive media. The am¬ 
plitude of this surface current is precisely that required 
to ensure that the rotation of the polarization of light 
reflected from this interface is zero, as is required by On- 
sager reciprocity in time-reversal invariant systems. Our 
result obtained from the Berry phase mechanism is the 
first example in which the surface current required by 
time-reversal symmetry emerges from a microscopically 
derived constitutive relation. 



FIG. 1. Momentum space representation of the Berry cur¬ 
vature mechanism for nonvanishing transverse current in a 
metal with inversion-symmetry breaking. The ellipsoid de¬ 
picts a typical Fermi surface with slices oriented perpendicu¬ 
lar to the optical wavevector. The Berry curvature, electron 
velocity and acceleration to first order in optical wavevector 
are illustrated for two points related by time-reversal. 


Tracelessness. As a preliminary step, we rederive the 
Berry-phase contribution to the gyrotropy tensor in a 
homogeneous time-reversal-symmetric material [18 . We 
use a ” Berry/Boltzmann” approach, in which the stan¬ 
dard calculation of linear response via the Boltzmann 
equation is augmented by the anomalous velocity term 
in Eq. [l] Solving the Boltzmann equation in the re¬ 
laxation time approximation yields the change in 

the distribution function to first order in the wave field 
E(r, t) = E exp (iujt — iq • r), 


( 3 ) 


where m 


5k = 


eE 


H(1/t — iu> + iq ■ v) ’ 


( 4 ) 


For simplicity, we concentrate in the following on the 
clean or high-frequency limit uor —>• oo. The semiclas- 
sical equations are valid as long as the frequency is be¬ 
low that of interband transitions and neglect electron- 
electron interactions (although incorporating electrons at 
the density-functional level is simple just as for the in¬ 
trinsic anomalous Hall effect). The current that arises 
from the anomalous velocity is given by[ 20 ], 


■I 


d 3 k 

(27T) 1 


f^Sk x St, 


where to first order in q, 


5k'. 


- eE(l + s _v)_ 


h \ co 

Substituting Eq. [ 6 ] into Eq. [5] we obtain 

- 2 r d 3 k 


j = ¥ Ex 


/(7^ + V) 




( 5 ) 


(6) 


( 7 ) 


The q-independent component of the integral in Eq. 
0 vanishes because time-reversal symmetry enforces 
ft( k) = — fi(—k). However, as Fig. 1 illustrates, the 
q-dependent term can be nonzero in the presence of IS 
breaking. An explicit example of a tight-binding Hamil¬ 
tonian with Berry curvatures of the required type was 
previously given [18]. The ellipsoid represents a typi¬ 
cal Fermi surface and the two parallel discs are slices of 
momentum space perpendicular to the wavevector of the 
light. Focusing on two representative points related by 
time-reversal, we see that the acceleration k to first or¬ 
der in q (shown as red arrow) is proportional to v and 
is therefore odd in k. Consequently the second term in 
the integrand of Eq. [7] is overall even and leads to a 
nonvanishing transverse current. 

Next, we re-express Eq. [7] in the standard form for the 
nonlocal constitutive relation, 


dE 

ji(u) = + 7 


(8) 


which relates the current to the first order of spatial 
derivative of the electric field US]. Using Eq. [7| 

lijk = JKjJ (2^y /(0) ^ Aufe ’ (9) 


where is the antisymmetric tensor. This response 
derived from the Berry curvature satisfies the condition 
7 iji = — 7 jn imposed by time-reversal symmetry [15, Tl . 

Because 7 ^ is antisymmetric, the gyrotropic response 
is usually expressed by its dual second-rank tensor g^, 
i.e., ji = —ieijkgidEjqi. Converting to this notation, 


e 


2 


9ij ~ ihu 



/ (1) = -V fe / (0) • (5k, 


( 10 ) 
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The trace of gij is given by, 


Y,9« 


ihuj (27t) 3</ ’ 


( 11 ) 


which is zero for the ground state SB or any other dis¬ 
tribution /depending only on energy. To see that 
integral over occupied states of ft • v vanishes even in 
the presence of monopole singularities in the Berry cur¬ 
vature, we write, fiv(k) = nds/dk ±, where n is normal to 
the surface of constant energy in momentum space and 
dk± is the separation between two such surfaces whose 
energy differs by ds. With this relation, the integral over 
occupied states can be written [22], 


calculations of a-quartz and trigonal Se [24], which iden¬ 
tify a contribution that is traceless within numerical er¬ 
ror (the cv part in that work’s notation). For a given 
material, measuring the trace of the gyrotropy tensor 
tests whether the Berry mechanism dominates other pos¬ 
sible contributions to the gyrotropic response, for exam¬ 
ple from the Bloch electron magnetic moment (spin [25] 
or orbital). 

The gyrotropic response that results from the Berry 
curvature is related to the question of the existence of a 
“chiral magnetic effect”, a current induced by a magnetic 
field in the presence of pairs of Weyl nodes (related to a 
triangle anomaly [23] l26ff30] ). According to the semiclas- 
sical theory [22] . the equilibrium current is 


/ (i^ /(0)v ' O = / dejjsn-h. (12) 



f^d 3 k 

(2tt) 3 


ft • V, 


(14) 


The integral is clearly zero in the absence of singularities 
in ft , as in this case V • ft = 0 for all k. However, 
the integral still is equal to zero [22: in the presence of 
singularities such as Weyl points, since 

/ d 3 k 

/ (0) V- ft = (/J.-£min)J2 qn ’ ( 13 ) 

which vanishes as the net monopole charge in the Bril- 
louin zone is zero because of lattice fermion doubling j23j . 

Tracelessness of gij has verifiable observable conse¬ 
quences: it is equivalent to the statement that the sum 
of the optical rotatory power measured along three prin¬ 
cipal axes is zero. This rule, derived on the basis of 
Berry-Boltzmann physics, goes beyond the constraints 
imposed by point group symmetry. There are fifteen 
crystal classes in which non-vanishing components of gij 
are allowed. Of these, eleven are chiral, indicating that 
all mirror symmetries are broken, and four have broken 
inversion symmetry but are not chiral. Point group sym¬ 
metry requires only these latter four classes to have trace¬ 
less gyrotropic tensors. Thus it would seem that for the 
other classes the observation of tracelessness would indi¬ 
cate the dominance of the Berry phase mechanism. 

There are hints that the Berry phase mechanism is ap¬ 
plicable to insulators as well as metals in the optical prop¬ 
erties of <a-quartz, one of the earliest and most studied of 
condensed matter chiral systems mmm- Point group 
symmetry applied to <a-quartz, which belongs to crystal 
class 32, requires only that (in the principal axis frame) 
two of the three diagonal elements of gij are equal, and 
the off-diagonal components are zero. Nevertheless, it is 
found experimentally that gn = #22 = = ( 1 / 2 )^ 33 , that 
is, the tensor is traceless over a broad frequency range 
that extends from visible to near-UV wavelengths. As 
it is extremely unlikely that this is accidental, there is 
evidence that at least in certain non-metallic systems a 
Berry phase-related mechanism is responsible for the gy¬ 
rotropic response. A hint is found in detailed ab initio 


and is therefore zero according to the argument presented 
above. However, the constitutive relation we have de¬ 
rived for the nonlocal current gives a closely related ex¬ 
pression for the current that accompanies a plane elec¬ 
tromagnetic wave. Consider a plane wave propagating 
along a principal axis of the crystal, which we take to be 
the z-direction. According to Eq. [lO] 

where we have used the Maxwell relation V x E = 
—<9B /dt. Thus the correct constitutive relation for the 
Weyl state is closely related to Eq. [14] but with the 
crucial differences that the response must be intrinsically 
anisotropic because of tracelessness of ftiVj, and the cur¬ 
rent must vanish if the magnetic field is static. 

Interfacial surface current. Combining Eq. [5] with 
Maxwell’s equations yields a difference in index of re¬ 
fraction for left and right circular polarizations, Sn± = 
n + — n_. The latter implies a rotation of the plane of lin¬ 
ear polarization with propagation through the medium, 
which is the phenomenon of NO A. At first glance, Sn± / 
0 would appear to predict polarization rotation on reflec¬ 
tion as well. The Fresnel formula for normal incidence 
reflection, yields a Kerr angle, 


However, the expectation that 5r± / 0 has been shown 
to violate the general reciprocity principle for electromag¬ 
netic fields interacting with time-reversal invariant media 
in equilibrium m- The seeming paradox is reconciled by 
an historically somewhat obscure [321134] (but recently 
rediscovered [35]) constraint on the nonlocal response 
functions imposed by time-reversal symmetry in media in 
which the gyrotropic coefficient varies in space, for exam¬ 
ple at an interface. While the validity of this constraint 
is not in question, there has yet to be a derivation of a 
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nonlocal constitutive relation in spatially varying chiral 
media that is consistent with time-reversal symmetry. 

A strength of the semiclassical approach is that com¬ 
bining real-space and momentum-space dependence is 
easier than with diagrammatic methods. We now include 
spatial variation of and show that the tracelessness 
of Qij is the crucial ingredient needed to obtain a fully 
consistent constitutive relation at the interface. 

Imagine that either by spatial variation of chemical 
composition, or some form of gating (see Fig. 2), a step 
in potential, V(z), is engineered that is sufficiently slow 
on the scale of the mean free path such that the semiclas¬ 
sical equations remain valid. We calculate the interfacial 
current in response to a plane electromagnetic wave with 
wavevector q z , to order E x V'(z). 

Including V(z), the equilibrium distribution function 
/becomes a function of z, leading to an extra term in 
the Boltzmann solution, = (df^/dz)5z, where, 


e (E x fi) 

h iuj 


z 


(17) 


is the distance travelled in the z-direction by a carrier in 
one cycle of the optical frequency. A detailed derivation 
of the component of to order V'(z)E is given in the 
Supplemental Material. The additional current arising 
from the spatial variation of /(°) is, 


j 



d 3 k df 
(27r) 3 dz 


(E X n) z v(k). 


(18) 


When the width of the interface is much less that the 
wavelength of the light, the relevant observable is the 
surface sheet current, 

B«W < 19 > 

where A/(°) is the change in /across the interface. 
Eq. ( fl9| ) corresponds to a constitutive relation for the 
interfacial current of the form, Ki = GijEj , where the 
surface conductance is given by, 

z 4 h, 

Gij = ifaj (2^ A/ efe ^ fe (k)^(k)- (2°) 

If the region z < 0 is emptied of carriers, such that we 
have an interface between gyrotropically active and inac¬ 
tive media, then A/(°) = /(°). The antisymmetric part 
of the surfuce conductance is 

1 p 2 r ^3 h, 

2 (G xy - G vx ) = — J —fW(Sl x v x + n y v v ), (21) 

or | ( G xy — G yx ) = g zz , by the tracelessness of gij. Fi¬ 
nally, we obtain for the antisymmetric part of the current 
response at the optically active/inactive interface, 

jx = gzz [d z + — S(fi)\Ey. ( 22 ) 


While the factor of 1/2 appearing in Eq. (22) can be 
shown to be required by time-reversal symmetry [21], it 
has not previously been derived from a microscopic or 
phenomenological theory. For uniaxial materials the con¬ 
stitutive relation Eq. [22|together with standard boundary 
conditions on the fields yields zero polarization rotation 
on reflection, as required by reciprocity. 



FIG. 2. Illustration of a slab of an acentric metal in 
which electrons are depleted underneath a gate electrode. In 
the presence of an electromagnetic wave, counterpropagating 
sheet currents appear at the interfaces betweeen optically ac¬ 
tive and inactive media. 


We note that in Eq. (22) time-reversal symmetry is 
preserved globally, but not locally. The conductivity ten¬ 
sor that describes the interfacial current violates Onsager 
reciprocity, as it is a local relation with antisymmetric 
off-diagonal components, z.e., G xy # G yx . Onsager reci¬ 
procity and time-reversal symmetry are restored only 
when considering the combined bulk and surface re¬ 
sponse. This behavior is reminiscent of 3D topological 
insulators, whose surface states have an odd number of 
Dirac fermions, which is impossible for a 2D time-reversal 
symmetric system in isolation. It would be worthwhile to 
understand possible additional electronic contributions 
to gyrotropy beyond the semiclassical static limit, as has 
been done for the chiral magnetic effect [361 in a Weyl 
semimetal model, where the “uniform” (not static) ef¬ 
fect is nonzero but not quantized. The Weyl semimetal 
TaAs @HZ], along with the similar candidate materials 
NbAs [37H39] and TaP [TO] 01], breaks inversion but its 
space group (Z4irad, #109) has point group 4mm, which 
does not allow optical activity. Either finding a differ¬ 
ent Weyl semimetal with lower symmetry, or lowering 
the symmetry of the TaAs family (e.g., by strain), would 
lead to a useful testbed for the Berry-phase contribution 
to gyrotropy as the magnitude of the Berry curvature is 
large near the Weyl points, leaving aside the possibility 
of open Fermi surfaces [42] . 

The results presented here are important for the emerg¬ 
ing field of topological photonics [43,. To date, this re¬ 
search has focused mainly on intricately fabricated meta¬ 
media in which response functions vary periodically on 
the scale of the optical wavelength. An example of a 
topological photonic state that can be created in this way 
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is an analogue of the quantum Hall effect [44], but inter¬ 
faces between conventional materials, which require less 
difficult fabrication, can also support topological inter¬ 
face states [45, [46] . However, in Ref. [45] the gyrotropic 
response is modeled as either a pseudoscalar, or a traceful 
tensor with a single diagonal component, both of which 
are excluded by our analysis. Thus, one implication of 
our findings is that future analysis of chiral/nonchiral in¬ 
terfaces should include the traceless property of gij. Fi¬ 
nally, an important open question is the relationship be¬ 
tween the interfacial photonic states generated by gij and 
the time-reversal protected electron conductance, g zz . 
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Supplemental Material for “Optical gyrotropy from axion electrodynamics in momentum space” 


For the specific example considered in the main text, a slowly varying scalar potential V(z) is added to the originally 
uniform material. The effects of V ( z ) come in two ways: the original distribution function will be inhomogeneous via 
a dependence on z, f°(z), and the semiclassical equations will be modified as 


r = v — k x $1, 

hk = —eE — er x B — V\z)z. 

With these we solve for / 1 which is the correction to /° in Boltzmann equation 

df . df , df 1 df 1 . df 1 ■ f 1 

o — Xi + ~^r + Y — x i + ~ -• 

OXj okj at oxj okj r 


(23) 

(24) 


(25) 


We will only consider terms linear in E field and keep up to the first order term in the spatial derivative. We also 
set r oo. 

First, consider the electric field component in x direction E x = Eq exp (iqz — iut). The solution 


f = 


df 


cE x 


+ 


df° eE x Q y 


dk x + iqv z ) dz ihuo 

can be seen directly by substituting it to the Boltzmann equation. We get 


df 1 . <9 2 / 0 


cE x 


dz 

sjl kz = _ 

dk z z dk x dk z h(—iuj + iqv z ) 


dk x dz + iqv z ) 

0 2 f e.E x 


df iqeE x 

V Z + ^7 - TZ-V Z , 

dk x —injjj 

(-V'(z)), 


d -L = + d Jl^ Vz - d JleEM 

dt hdk x dk x ihuj hdz y ' 

df df 

~d^ i= Wz^ Vz+eEx ^ v) ' 

df. f0 

dk x kx = hdk~f~ eEx + eVzBy) ’ 
df . df 


The summation of these terms is zero. 
Similarly, for E y and E z the solutions are 


f = 


df 


eE v 


+ 


df 


0 ~SEyQ X 


dky h{—iuj + iqv z ) dz iftio 


and 


f = 


df 


eE, 


dk z h(—iu> + iqv z ) 


(26) 


(27) 


(28) 


(29) 


respectively. 

With these f 1 , we compute the coefficients 7 and A in the standard expression of current j for optical rotation: 

ji = ^ijEj T 72 j 1 V 1 Ej T Ej V [Xiji. (30) 

There is a minor subtlety in computing the current at finite r (one that also appears in the semiclassical treatment 
of the anomalous Hall effect) that we discuss in the appendix. Note that \ji can be divided into two parts: the 
symmetric part A s and the antisymmetric part X an . The symmetric part V^Af^ can be reorganized into e^-. 
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The first term in / 1 contributes to 7 because the steaming of electrons into or out of a certain k space region will 
change the polarization. The second term contributes to A through the normal effect of a shift of the Fermi surface. 
Note that, at the order we are concerned with, we should also include effects from the magnetic field but those turn 
out to cancel out via the tracelessness relation. 


Below we list all the components explicitly. 

r )yxz 
7 xyz — 

7 xzz = 

7 zxz 

7 zyz = 
lyzz = 


— 1 

r d 3 k 

itvjj J 

(27r) 3 

-e 2 

f d 3 k 

iftoj J 

(27r) 3 

— 1 

f d 3 k 

ihhj J 

(2n) 3 

—e 2 

f d 3 k 

ihuj J 

( 2 t r ) 3 

— 1 

f d 3 k 

ihuo J 

(27r) 3 

-e 2 

f d 3 k 

ihhj J 

(27r) 3 


foV z Qz 1 

foV z Qz 5 
foVz^ly, 


X yxz 


^xyz — 


^xzz — 0 , 


•fc. . / /o_\3 foVz^yt ^zxz — . fc . . 


^zyz — 


foV z ^x •> 

foV z ^x, \zz *= 0. 


—e 2 

r d 3 k 

i/kj J 

( 27r) 3 

_ 1 


ihuj J 

( 27 t ) 3 

0, 

—e 2 

r d 3 k 

i/kj J 

( 27r) 3 

— 1 

f d 3 /c 

ihuj J 

( 27 t ) 3 


f0Vy£ly, 

foV x ^x •> 

foV z Qy 5 

foV z ^x •> 


(31) 


We should have other components of A like \ xxz but these terms clearly belong to A s . 

The components of 7 above have the desired symmetry properties 7^7 = — 7 ju and it can also be shown that 
A?jl = First, by using the tracelessness relation 


Kxz - 2 (V Kvz) - 2 iHu J (2?r) 3 foV z n 


is half of 7 yxz . The remaining relations can be seen straightforwardly, like 

d 3 k 


\an = I - e2 

2 ihu 


/ 


— / 773/0^ 


(32) 


(33) 


which is half of 7 zxz . 



